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ABSTRACT 

We report the discovery and characterisation of a new M-dwarf binary, with 
component masses and radii of Mi = 0.244to;oo3= O-^Sllo pog ^ 0 , M 2 = 
0.179to'ooi-^0) ^2 = 0.218to;oii-^0) orbital period of ~ 4.1 days. The M-dwarf 
binary HATS551-027 (LP 837-20) was identified as an eclipsing binary by the HAT¬ 
South survey, and characterised by a series of high precision photometric observations 
of the eclipse events, and spectroscopic determinations of the atmospheric parame¬ 
ters and radial velocity orbits. HATS551-027 is one of few systems with both stellar 
components lying in the fully-convective regime of very low mass stars, and can serve 
as a test for stellar interior models. The radius of HATS551-027A is consistent with 
models to Icr, whilst HATS551-027B is inflated by 9% at 2cr significance. We measure 
the effective temperatures for the two stellar components to be Teff.i = 3190 ± lOOK 
and Teff ,2 = 2990± IlOK, both are slightly cooler than theoretical models predict, but 
consistent with other M-dwarfs of similar masses that have previously been studied. 
We also measure significant Ha emission from both components of the binary sys¬ 
tem, and discuss this in the context of the correlation between stellar activity and the 
discrepancies between the observed and model temperatures. 
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data gathered with the 6.5 meter Magellan Telescopes located at 
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1 INTRODUCTION 

Double-lined M-dwarf eclipsing binaries are natural labora¬ 
tories for stellar astrophysics, providing the most precise, 
model independent, measurements for the fundamental stel¬ 
lar parameters of low mass stars. They have been used ex¬ 
tensively to test our theoretical understanding of stellar in - 
teriors (e.g. lChabrier k, Baraffell 19951 : [Torres Ril^l2002l 'l. 
Very low mass stars (VLMS), with masses below 0.35 Mq, 
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are thought to have fully convect ive, adiabatic interiors 
fe.g. IChabrier fc Baraffelll997l . [^0di . Models have been un¬ 
able to match the observations of VLMS, with measured 
temp eratures cooler and radii larger than models predict 
(e.g. Torres fc Ribad 200^: Rib^ 200d: Torres et al'l 2010l : 


iFeiden fc Chaboveill2ni2l : ISnada et al.l lTorredl2ni3h . 

Double-lined M-dwarf binaries in the VLMS mass 
regime, where the stellar parameters have been ac¬ 
curately determined, are rare. The CM Draconis sys¬ 
tem is the most well studied, with masse s and radii 
determined to bet t er than 1% precis ion (iLac 3 1 19771: 
iMetcalfe et al.iri99d : [Morales et ahlfiood l. The radius dis¬ 
crepancy between models and observations have been 
well documented fo r CM Draconis at the 2-10% 
level dT orres fc R. bad l2002l: iMacDon ald fc MullanI l2012l : 


ISpada et al.l 20131 : Feiden fc Chabov^ 20l4 l. depending on 
the treatment of metallicity, mixing length, and mag¬ 
netic inhibition in the models. More recently, the systems 
KOI-126, Kepler-16, LSPM J1112-b7626 and WTS 19g- 
4-02069_( [Carto et al.l 1201 ll : IPovle et al.l l 201 ll : |i rwin et ^ 
I2OIII : iNefs et ^ 2OI3I I have been reported, with masses and 
radii of the low mass stellar components measured to bet¬ 
ter than 2% precision. Of these systems, only the radii of 
KOT126B and C were found to be in agreement with the 
models of slightly super-solar metallicity, wh ile the other 
VLMSs are inflated compared t o the models llFeiden et al.l 
201l|^ Spa da fc Demaraua I2OI2I : iFeiden fc Chaboverl 20 14 


Nefs et ^ 2013l~). 


A number of ideas have been put forward to explain the 
radius and temperature discrepancies between observations 
and the models. Increased metallicity and missing opacity 
may account for the larger measured ra dii of some systems 
llBerger et al.l[2003 : iBurrows et al.ll201ll l. In particular, the 
KOI-126 B and C companions can be well modelled by as - 
suming a metallicity of [Fe/H1 = -|-0.15 (|Feiden et al.l[20lll l. 
However, iLopez-MoralesI (I 2 OO 7 I I found that a metallicity - 
radi us relationship exists for si ngle stars, but not for bina¬ 
ries. ISpada fc Demarauel (I2OI2I I suggested increased metal¬ 
licity cannot fully ac count for th e mod el discrepancy for the 
CM Draconis binary. IZhou et al.l (120141') found no correlation 
between metallicity and the measured radii of the VLMS 
sample, suggesting that the effect of metallicity variation on 
radius is smaller than the present measurement uncertain¬ 
ties. 

The spin-up of M-dwarfs in binary systems, lead¬ 
ing to higher magnetic activities, ha s also been put for¬ 
ward to resolve the d i screpancies (e.g . jLo pez-MoralesI 120071 : 
IChabrier et al.llioOTl l. IChabrier et al.f (l2007l l suggested that 
strong magnetic fields can inhibit convection and heat 
flow within the stars; increased spot coverage from higher 
magnetic activities will also modify the evolution of ac¬ 
tive low mass stars, leading to different radii than in¬ 
active stellar models. This is supported by the empiri¬ 
cal correlation between X-ray and Ha activity leve l s and 
M-dwarf radii (iLopez- MoralesI l2007l : [Morales et al.l I2OO8I : 


Stassun_et_^ 201^. M odels of the CM Draconis system by 


MacDonald fc MullanI (l2012l . l2014h found that magnetic in¬ 


hibition of convection can explain the radius discrepancy. 
However. IFeiden fc Chaboverl J2014h modelled the Kepler-16 
and CM Draconis systems, and found it difficult to sustain 
the strong magnetic fields required to inhibit convection and 
inflate the stars to the observed radii. 


Recently, accurate M-dwarf parameters have also be¬ 
come important for the characterisation of the planets they 
host. The larger mass, radius, and luminosity ratios between 
M-dwarf hosts and planets make them enticing targets for 
discovery and follow-up characterisation observations. The 
discrepancies in theoretical modelling of M-dwarfs resulted 
in a need for empirical relationships, derived from well char¬ 
acterised M-dwarf measurement s, in order to characterise 
the star and planet s ystems (e.g. [Johnson et al.l[201ll . l2012l : 
iHartman et al.ll2014l ). 

In this paper, we present the discovery and charac¬ 
terisation of a new eclipsing, double-lined, M-dwarf bi¬ 
nary, HATS551-027 (LP 837-20), with both stellar com¬ 
ponents firmly within the fully convective mass regime of 
VLMSs. HATS551-02 7 has pr eviousl y been identifie d as a 
high proper motion (iLuvted 1 19791 : fsalim fc Gouldl 120031) 
M4.5 dwarf jReid et al.ll2003l) from astrometric and low res¬ 
olution spectroscopic surveys of NLTT high proper motion 
stars. As part of our eff ort to characterise the V LMS popula- 
tion JZhou et al.ll20lj ). the HATSouth survey jBakos et al.l 
I2OI3I) identified and characterised HATS551-027 as a low 
mass eclipsing binary. We derive precise fundamental prop¬ 
erties of the system via the discovery light curves and a series 
of photometric eclipse follow-up observations and spectro¬ 
scopic orbit measurements, described in Sections [2] and [S] 
The HATS551-027 system is then discussed in the context 
of other VLMS binaries in Section U] 


2 OBSERVATIONS AND ANALYSIS 
2.1 HATSouth Photometric Detection 


The eclipses of HATS551-027 were first identified by ob - 
servations from the HATSouth survey dBakos et al.ll2oi3 ). 
HATSouth is a global network of identical, fully robotic tele¬ 
scopes, providing continuous monitoring of selected 128 deg^ 
fields of the southern sky. A total of 16622 observations 
of HATS551-027 were obtained from HATSouth units HS- 
1, HS-2 in Chile, HS-3, HS-4 in Namibia, and HS-6 in 
Australia from September 2009 to September 2010. Each 
HATSouth unit includes four 0.18 m f/2.8 Takahasi astro- 
graphs each with a Apogee 4Kx4K U16M Alta CCD cam¬ 
era. The images have a field of view of 4 x 4°, with a pixel 
scale of 3.7” pixel”^. The observations are performed in the 
Sloan-r' band, at 4 minute cadence. HATSouth photome¬ 
try is d etrended using the External Parameter Decorrelation 


fEPD. iBakos et ^ 2007) and Trend Filtering Algorithm 


(TFA, iKovacs et al. 1^0^ techniques. The transit candi¬ 
dates seajch_Js_pCTfomi with the Box-fitting Least Squares 
(BLS li^vacs et 3 II 2 OO 2 I ) analysis. Details of the HATSouth 
reduct ion a nd analysis pr o cesses can be found in iBakos et al.l 
1I2OI3I) and IPenev et al.l 1I2OI3I) . The HATSouth discovery 
light curve for HATS551-027 is displayed in Figure [T] the 
photometric data can be found in Table [5] 


2.2 Photometric Follow-up 

We performed photometric followup of the primary and sec¬ 
ondary eclipses of HATS551-027. The observations are de¬ 
scribed below, summarised in Table [T] with the light curves 
plotted in Figure [2] and found in Table [S] The specific set 
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Figure 1. HATSouth discovery light curves of HATS551-027. The observations are plotted in grey, best fit model from Section 13.11 
plotted in red. Close-ups of the primary and secondary eclipses are plotted in the lower panels. 
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Figure 2. Follow-up photometric observations for the primary (left panel) and secondary eclipses (right panel) of HATS551-027. The 
observations are plotted in black, best fit model from Section l3. 1 1 plotted in red. The uncertainties shown have been scaled up to force a 
reduced = 1 where necessary. 
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Table 1. Summary of photometric observations 


Facility 

Dates 

Number of Images 

Exposure Time (s) 

Filter 

HATSouth Network 

2009 Sep - 2010 Sep 

16622 

240 

r' 

FTS / Merope 

2012 Dec 12 

124 

60 

i' 

Swope / SITeS 

2013 Feb 26 

78 

30 

i' 

FTS / Merope 

2013 Mar 20 

114 

60 

i' 


Table 2. HATS551-027 differential photometry 


BID 

Flux 

A Flux 

Instrument 

Filter 

2455083.75892 

0.9999 

0.0236 

HATSouth“ 

r' 

2455083.76216 

1.0084 

0.0192 

HATSouth 

r' 

2455083.76552 

1.0514 

0.018 

HATSouth 

r' 

2455083.76873 

0.9989 

0.0209 

HATSouth 

r' 

2455083.7721 

0.9632 

0.0221 

HATSouth 

r' 

2455083.77532 

0.9487 

0.0204 

HATSouth 

r' 

2455083.77868 

1.0744 

0.0173 

HATSouth 

r' 

2455083.78191 

0.9683 

0.02 

HATSouth 

r' 

2455083.78528 

1.0118 

0.0196 

HATSouth 

r' 

2455083.78849 

1.014 

0.0188 

HATSouth 

r' 


Note. - This table is available in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form 
and content. 

“For the HATSouth light curve these magnitudes have been detrended using the EPD and TFA procedures prior to fitting a transit 
model to the light curve. Primarily as a result of this detrending, but also due to blending from neighbors, the apparent HATSouth 
transit depth is ~ 90% that of the true depth in the Sloan r filter. 


Table 3. Summary of reference stars for follow-up photometry 


Reference star 

RA 

DEC 

1 mag 

FTS 2012-12-12 / 2013-03-20 



1 

-1-05:44:56.682 

-24:57:43.844 

12.1 

2 

-1-05:44:48.663 

-24:56:55.899 

14.9 

3 

-1-05:45:03.111 

-24:56:18.213 

15.0 

4 

-1-05:44:59.082 

-24:54:44.068 

15.7 

5 

-1-05:45:02.925 

-24:54:28.156 

15.8 

6 

-1-05:44:56.505 

-24:54:12.312 

16.3 

7 

-1-05:44:49.306 

-24:56:45.252 

16.2 

8 

-1-05:45:06.857 

-24:55:22.837 

16.6 

9 

-1-05:45:03.468 

-24:56:38.683 

17.4 

10 

-1-05:44:51.649 

-24:54:24.864 

17.0 

11 

-1-05:45:02.931 

-24:55:32.118 

15.3 

Swope 2013-02- 

■26 



1 

-1-05:45:29.156 

-24:49:56.572 

12.7 

2 

-1-05:45:23.524 

-24:51:07.124 

13.0 

3 

-1-05:44:32.656 

-24:54:16.506 

13.4 

4 

-1-05:44:30.314 

-24:54:57.665 

13.5 

5 

-1-05:45:10.934 

-25:02:41.212 

13.2 


of reference stars and their magnitudes for differential pho¬ 
tometry are listed in Table O 

2.2.1 Faulkes Telescope South 2m / Merope 

Two secondary eclipses of HATS551-027 were observed by 
the Merope camera on 2 m Faulkes Telescope South (FTS), 
at Siding Spring Observatory, on 2012 December 12 and 


2013 March 20. FTS/Merope has a field of view of 4.7 x 4.7', 
and a pixel scale of 0.139” pixel~'^, read-out using 2x2 bin¬ 
ning mode. The observations were obtained in the i-band, 
with exposure times of 60 seconds. The telescope was defo- 
cused to reduce the effect of intra- and inter-pixel variations 
on the final light curve, and to avoid saturation of target and 
reference stars. A total of 124 observations were obtained on 
2012 December 12, and 114 on 2013 March 20. 

Bias and flat field reduced images, automatically gen¬ 
erated by the LCOGT reduction pipeline, were used for the 
photometric extraction. The raw light curves for the target 
and reference stars were extracted using Source Extractor 
llBertin fc Arnout^ll996h . via multiple fixed circular aper¬ 
tures, and after an interpolated background map was sub¬ 
tracted. Eleven reference stars of magnitudes / = 12.1-17.4 
were used for the relative photometry. The best aperture was 
chosen to minimise the out-of-transit light curve scatter. A 
second order polynomial was fitted to the out-of-transit re¬ 
gions of the light curves to remove small, long-duration vari¬ 
ations, likely due to changing airmass, before model fitting 
was performed. 

2.2.2 Swope Im / SITeS 

A near-complete primary eclipse of HATS551-027 was ob¬ 
served by the SITe#3 camera on the Swope 1 m telescope 
at Las Campanas Observatory, Chile, on 2013 February 26. 
Swope/SITe3 has a field of view of 14.8 x 22.8', a pixel scale 
of 0.435” pixel”^, and detector size of 2048x3150 pixels. The 
observations were performed in the f-band, with 30s expo¬ 
sures. A total of 78 exposures were obtained. Photometric 
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extraction were performed on the reduced frames using the 
FITSH package (H a IH). Five reference stars of magni¬ 
tudes I = 12.7-13.5 were used for the relative photometry. 
We note that no pre-ingress baseline was recorded, which 
may have an adverse effect on the derived timings and asso¬ 
ciated uncertainties. 


2.3 Spectroscopic Follow-up 

We performed a series of spectroscopic follow-up observa¬ 
tions at low, medium, and high resolution, to determine the 
stellar properties and orbit of the system. The observations 
are summarised in Table[4l and the radial velocities are plot¬ 
ted in Figure [3] and presented in Table [S] 


2.3.1 ANU 2.3m / WiFeS 


Spectroscopic observations at low and medium resolutions 
were performed us ing the Wide Field Spectrograph (WiFeS, 
iDopita et al.ll20(^ ) on the ANU 2.3m telescope, located at 
Siding Spring Observatory, Australia. WiFeS is an image 
sheer integral field spectrograph, with slitlets of 1” width 
in spatial coverage. The WiFeS observations provided the 
spectroscopic classification and radial velocities (to the ~ 
2kms"i level) of HATS551-027. 

For spectral classifications, an observation at low reso¬ 
lution of A/A A = R = 3000 was obtained using the B3000 
and R3000 gratings on the blue and red arms of the spectro¬ 
graph, with the RT560 dichroic, providing a flux calibrated 
spectrum of the object over 3500-9000 A. For the R = 3000 
observations, the stellar flux is summed over the three slitlets 
together to produce the final spectrum . Flux calibrations are 
performed according tolB esselll lll999ll. using flrrx standard 
stars from iHamuv et aLl i 1992fl and lBesselll (Il999l l taken on 
the same night. The observational techniques, data reduc¬ 
tion and flux calibration processes for t he WiFeS low resolu - 
tion observations are fully described in iBavliss et al.l (l2013ll . 

To measure the radial velocity orbit, we then ob¬ 
tained 16 medium resolution {R — 7000) observations from 
September 2012 to August 2014. The observations were per¬ 
formed using the R7000 grating and the RT480 dichroic, 
giving a velocity resolution of 21.6 kms”*^ pixeR*^ over the 
wavelength range 5200-7000 A, corresponding to the photo¬ 
metric R band. The wavelength calibrations were obtained 
from bracketing Fe-Ne-Ar arc lamp exposures. In addition, a 
first order correction to the wavelength calibration is applied 
using the telluric Oxygen B lines over 6882-6906 A. For ra¬ 
dial velocity measurements, we reduce the spectra from each 
WiFeS slitlet separately, such that three reduced spectra 
from three slitlets are produced per exposure. The reduc - 
tion process is described in detail in iBavliss et al.l (l2013ll . 
We adopted the Ha-derived velocities for analysis in Sec¬ 
tion [33] 

Radial velocities of both components of the binary were 
first measured by fitting the Ha emission feature, which ex¬ 
hibits clear and distinct signals from both stellar compan¬ 
ions. For each spectrum, we fit the Ha feature with two 
Gaussian components. The best fit values are obtained via a 
down-hill simplex minimisation using the Python fmin func¬ 
tion in scipy package, and the errors explored via an MCMC 
analysis using the emcee affine invariant ensemble sampler 


llForeman-Mackev et al.l[2013l l. The final velocity error from 
each exposure is a quadrature combination of the mean cen¬ 
troid errors from the MCMC fit and the scatter in velocity 
from three slitlets. The Ha emission features and Gaussian 
fits for each observation are presented in Figurejd] The obser¬ 
vations from each night are plotted by three black solid lines, 
for the spectra from the three slitlets. The best fit Gaussians 
for each spectrum are also plotted by the dotted line, in red 
for the primary component, blue for the secondary compo¬ 
nent. The velocities from WiFeS are presented in Figure [3] 
and Table (5) 

We also tried cross correlations between the target spec¬ 
tra and an M4V template from a WiFeS observation of 
the standard star LP 816-60. The radial velocity standard 
cross correlations provided blended cross correlation func¬ 
tions (CCFs), which did not resolve the two spectral com¬ 
ponents as well as the Ha feature. We derived radial ve¬ 
locities from selected WiFeS R7000 observations taken at 
phase quadrature by simultaneously fitting a double Gaus¬ 
sian to the CGF profile, presented in Figure [3] and Table (5) 
We also tried cross correlating against spectral template s 
from the BT-Settl model atmospheres (I Allard et al.|[201^ '). 
but we could not resolve the two spectral components. 

The WiFeS cross correlations were also used to measure 
light ratio between the two stellar components (L 2 /T 1 ). We 
fit a double Gaussian to the GCFs from the spectra, and 
derived a light ratio of 0.53 ± 0.06 from the R7000 grat¬ 
ing, equivalent of the photometric R band. We also obtained 
three additional exposure over the wavelength range 6800- 
9100 A at i? = 7000, using 17000 grating and the RT615 
dichroic, corresponding to the photometric I band. This ex¬ 
posure was used to measure the light ratio between the two 
components at the wavelength of the follow-up photomet¬ 
ric observations. The 17000 observation gave a light ratio of 
0.50 ± 0.01, consistent to the R band light ratio. 

2.3.2 MPG 2.2m / FEROS 

Four observations of HATS551-027 were obtained using the 
fibre-fed FEROS spectrograph on the MPG 2.2m telescope 
at La Silla Observatory, Ghile , over the nights 2013-03-24 
to 2013-03-27. FEROS provides a resolution of R = 48000, 
descriptions of the FEROS observations and r eduction pro¬ 
cedures can be found in ijPenev et al.l I2OI3I : [Jordan et al.l 
I2OI4I : iBrahm et al.l 120151 1 . The observations are cross cor¬ 
related with an M-dwarf template. The two components of 
the binary are fully separated in the CGF, and the veloc¬ 
ities are presented in Figure [3] and Table [3] We also mea¬ 
sure the Gaussian peak heights of the CCFs to estimate the 
light ratio for the two components of the system. Because of 
the faintness of the target star, we find that two orders in 
the wavelength range 6475-6665 A reliably showed the two 
CCF peaks in all of the exposures. The average and stan¬ 
dard deviation CCF peak heights over the 4 exposures was 
0.42 ±0.07. 

2.3.3 Magellan 6.5m / PFS 

Four consecutive observations of HATS551-027 were ob- 
tained using the C arnegie Planet Finder Spectrograph (PFS, 
ICrane et al. l201Clh on Magellan H, Las Campanas Observa¬ 
tory, Chile, on the night of 2013-02-21 UT. The observations 
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Table 4. Summary of spectroscopic observations 


Facility 

Date Range 

Number of Observations 

Resolution 

S/N“ 

ANU 2.3m / WiFeS 

2012 Sep 8 

1 

3000 

95 

ANU 2.3m / WiFeS 

2012 Sep 8 - 2014 Aug 14 

16 

7000 {R band) 

20 

MPG/ESO 2.2m / FEROS 

2013 Mar 24 - 27 

4 

48000 

35 

Magellan 6.5m / PFS 

2013 Feb 21 

4 

76000 

90 

ANU 2.3m / WiFeS 

2014 Sep 6 - 2014 Oct 10 

3 

7000 (/ Band) 

60 


“ Average signal-to-noise ratio per resolution element, at approximate peak. 



Phase 

Figure 3. The radial velocity orbit of HATS551-027A (red) and B (black). The best fit models from Section |3T] are plotted by the red 
and black curves for components A and B respectively. 


were taken through the iodine cell, with a 0.5” x 2.5” slit 
and using 2x2 binning in slow read-out mode. 

We use the PFS spectrum to derive the single epoch 
relative radial velocities and light ratio of the two stellar 
components. The spectra from the four observations were 
averaged, and cross correlated against a Tefj = 3100 K, 
log q = 5.0 solar m etallicity BT-Settl synthetic spectrum 
ll Allard et al.ll201^ . We used five orders outside of the Io¬ 
dine absorption region (6114-6574 A) for the cross correla¬ 
tion. The two stellar components are resolved and fully sep¬ 
arated in the CCF. We simultaneously fit two Gaussians to 
the CCF, deriving the light ratios and the velocity centroids. 
The velocities from the five orders are weight averaged to 
arrive at the final relative velocity between the two stellar 
components. The light ratio I/ 2 /T 1 is 0.51 ±0.04, with the 
error as the standard deviation between the five orders. We 
also tried 3000 K and 3200 K synthetic templates, yielding 
the same light ratio to within ±0.01. The light ratio mea¬ 
sured from PFS is also consistent to ~ Icr of that measured 
from the FEROS spectra and WiFeS spectra. We adopt the 
light ratio measured from the high signal-to-noise PFS ob¬ 
servations for our modelling in Section r3.ll 


3 SYSTEM PARAMETERS 

3.1 Global modelling of light curves and radial 
velocities 

We perform a global modelling of the HATSouth discovery 
light curves, follow-up photometry, and spectroscopic radial 
velocities to derive the system parameter s for HATS551- 
027. W e model the light curves according to lNelson fc Da^^ 
( 197^, using a modified version of the JKTEBOP code 
( Popper fc EtzellIl98ll : ISouthworth et ahl |2004|) . The radial 
velocities are modelled using Keplerian orbits. 

The free parameters of the fit are period P, reference 
transit time to, radius ratio R 2 /R 1 , normalised radius sum 
(i?i±7?2)/a, line-of-sight inclination i, light ratio L 2 /L 1 , or¬ 
bital eccentricity parameters ecoscu and esinw, and orbital 
semi-amplitudes Ki, K 2 for the two comp onents. We assign 
a Gaussian prior to I/ 2 /±i (following e.g. lirwin et al.ll201lh 
from the PFS light ratio of 0.51 ± 0.04 from Section 12.3.31 
Without the prior, the radius ratio becomes degenerate with 
the light ratio, and is difficult to constrain. We assume uni¬ 
form priors for all other free parameters. The quadratic 
limb dark ening coefficien ts for both components are fixed to 
that from IClare 3 (l2000l) using the PHOENIX models, and 






















































An M-dwarf binary from HATSouth 7 


Table 5. HATS551-027 Radial Velocities 


HJD 

RVi (kms-i) “ 

(tRVi (kms 

RV 2 (kms-l) 

0 -RV 2 (kms-i) 

Instrument 

2456179.31471 

-0.8 

6.9 

103.7 

11.3 

WiFeS (Ha) 

2456179.31471 

-25.8 

3.3 

70.9 

6.7 

WiFeS (CCF) 

2456180.28765 

41.3 

10.5 

35.5 

14.6 

WiFeS (Ho) 

2456198.23861 

79.5 

12.1 

32.0 

19.7 

WiFeS (Ho) 

2456198.23861 

55.1 

2.5 

-9.9 

4.7 

WiFeS (CCF) 

2456231.09138 

44.7 

13.5 

32.3 

17.4 

WiFeS (Ho) 

2456375.98307 

50.9 

9.3 

28.9 

14.2 

WiFeS (Ho) 

2456378.02998 

51.3 

10.4 

31.8 

16.2 

WiFeS (Ho) 

2456378.97110 

-1.6 

3.1 

100.4 

4.6 

WiFeS (Ho) 

2456378.97110 

-18.9 

2.1 

74.3 

4.3 

WiFeS (CCF) 

2456380.99658 

84.5 

6.0 

-18.4 

11.4 

WiFeS (Ho) 

2456380.99658 

66.2 

2.1 

-26.8 

6.2 

WiFeS (CCF) 

2456876.30644 

2.4 

3.1 

102.0 

6.6 

WiFeS (Ho) 

2456876.30644 

-18.3 

2.8 

75.3 

5.1 

WiFeS (CCF) 

2456877.30697 

46.9 

11.3 

40.0 

16.7 

WiFeS (Ho) 

2456878.31132 

85.1 

7.3 

-9.8 

12.1 

WiFeS (Ho) 

2456878.31132 

63.3 

4.8 

-26.8 

6.2 

WiFeS (CCF) 

2456879.30918 

48.7 

14.2 

38.1 

18.7 

WiFeS (Ho) 

2456880.30174 

7.1 

2.9 

101.2 

5.5 

WiFeS (Ho) 

2456880.30174 

-15.8 

1.7 

76.6 

4.7 

WiFeS (CCF) 

2456881.28823 

50.1 

14.9 

48.1 

22.5 

WiFeS (Ho) 

2456883.25386 

56.3 

16.4 

19.8 

24.8 

WiFeS (Ho) 

2456883.25386 

50.0 

1.6 

-0.5 

5.6 

WiFeS (CCF) 

2456884.25377 

9.5 

4.5 

95.6 

6.7 

WiFeS (Ho) 

2456884.25377 

-8.5 

1.6 

71.5 

3.8 

WiFeS (CCF) 

2456375.53799 

-5.8 

0.2 

65.7 

0.6 

FEROS 

2456376.54952 

53.8 

0.2 

-16.1 

3.7 

FEROS 

2456377.53960 

56.6 

0.3 

-20.1 

4.8 

FEROS 

2456378.52904 

-1.5 

0.3 

60.4 

0.5 

FEROS 

2456344.60054 

3.2 

0.3 

-93.8 

0.9 

PFS 


“ The zero-point of these velocities is arbitrary. An offset is added to each instrument independently in the global fitting. 


given in Table (S] At each iteration, the modelled mass ratio 
(g = K 1 IK 2 ) is used to calculate the ellipsoidal variabilit y 
expected in the light curve. Following [Morales et al.l ll2009ll . 
we assign fixed reflection albedos for both components of 
0.5, and gravity darkening coefficients of 0.2. The instru¬ 
ment zero points for each dataset are not included as free 
parameters in the global model, but rather fitted for at each 
iteration independently. We often find the HATSouth dis¬ 
covery fight curves to be diluted due to the TFA detrending 
that was applied, as well as third-light blending from nearby 
neighbours due to the larger pixel scales of the instrument. 
The model transit shape is significantly constrained by the 
follow-up observations, and can be used to correct for the 
dilution of the HATSouth fight curves. At each iteration, 
we adjust the transit depth of the HATSouth light curves 
to t he tested model b y fitting for a third light component 
(e.g. IZhou et ^ lE^. We adopt the WiFeS radial veloci¬ 
ties derived from the Ha emission line in our analysis, since 
the two stellar components are better resolved in Ha than 
in the CCFs. 

The best fit parameters and the associated uncertain¬ 
ties are explored via a Markov chain Monte Carlo (MCMC) 
analysis, using emcee. For each MCMC run, we initialise 100 
walkers over 2500 steps each. The walkers are initialised with 
randomised starting positions around their respective ex¬ 
pected best fit parameter values. Examination of the MCMC 


chain shows that convergence is generally reached after 500 
- 1000 steps, depending on the initial conditions. We then 
re-run the MCMC with walkers initialised around the best 
fit parameters from the first run. We inflate the per-point 
measurement uncertainties such that the reduced for each 
dataset is at unity with respect to the best fit parameters 
derived from the first MCMC run. This allows other sources 
of error unaccounted for in data reduction to be included 
in the MCMC analysis. We found that only the FTS ob¬ 
servations required their measurement uncertainties to be 
inflate, by 2.4x and 2.lx for the 2012-12-12 and 2013-03-20 
observations respectively. 

3.1.1 Accounting for Spots 

Late M-dwarfs are known to be active, and the presence of 
chromospheric spectral emission features of HATS5 51-027 
suggest that it is n o exception. A Lomb-Scargle ILS. iLoinbI 
Il97d : IScarglelll982lf analysis of the out-of-transit HATSouth 
light curve residuals from Section 13.11 shows a highest peak 
at 4.24 days, with peak FWHM of 0.08 days, compared to 
an orbital period of 4.08 days. Although when we break the 
light curve into four segments of 4000 points each, we cannot 
demonstrate the consistency of the peak location. 

Star spot activity affect the observations and model re¬ 
sults in several ways. First, they induce periodically modu- 
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Wavelength (A) 


Figure 4. The Ha emission feature from the WiFeS R = 7000 
observations. For each exposure, the spectra from the three slitlets 
are plotted individually, in black. The best fit profiles are plotted 
in red for the primary component, and blue for the secondary 
component. Each exposure is arbitrarily offset in flux by 2 for 
clarity. 


lated out-of-transit variations in the HATSouth light curves. 
We correct this by modelling the cosine filtered light curve 
in the MCMC analysis. We model the modulation in the 
HATSouth light curve residuals using a set of linearly com¬ 
bined cosine functions, following Equation 1 of iHuang et al.l 
(|2013h . The HATSouth light curve residuals, with transits 
and eclipses masked, are smoothed by a 50 point-width box¬ 


car, and fitted with the cosine filtering model with minimum 
period of P/2 and maximum period of twice the observa¬ 
tional baseline. The light curve residuals and cosine filter¬ 
ing model for a segment of the HATSouth observations are 
shown in Figured The cosine filtering model has a standard 
deviation of 0.006, and peak-to-peak amplitudes of ~ 2%. 
For comparison, the modelled out-of-transit variations orig¬ 
inating from reflection and ellipsoidal effects is at the 10~^ 
level. 

Variations in spot coverage also alters the total flux 
of the system, leading to variations in the eclipse depth. 
This affects the parameters derived from the high precision 
follow-up photometry, such as R 2 /R 1 and {R 2 +R 1 ) /a. Since 
the follow-up observations were performed 800 days af¬ 
ter the end of the HATSouth continuous observations, they 
cannot be corrected by emp loying spot model f its to the 
out-of-transit variability fe.g. [Morales et aHl2009l ). Instead, 
we create 100 different sets of the follow-up light curves, 
with each light curve scaled by a random factor drawn from 
a Gaussian distribution with ^ = 1.0 and cr = 0.006 (the 
standard deviation of the cosine filtering model). We then 
perform a global analysis detailed in Section 13.11 using the 
scaled sets of the follow-up light curves. The MCMC chains 
from each of the 100 separate analyses are combined to pro¬ 
vide the final system uncertainties. 

In-transit star spot crossing events modify the shape of 
the transit light curve, and can be identified from the corre¬ 
lated red-noise in the fit residuals. Whilst we cannot correct 
for the spot crossing events, we can account for their effects 
in the uncertainties of the model fit. To account for corre¬ 
lated noise, w e treat the follow-up photometry as a Gaussian 
process (e.g. ICibson et al.ll2012 l. We form the covariance 
matrix E, of dimension N x N, where N is the number of 
data points per follow-up light curve. The elements of the 
covariance matrix is defined as 

T,ij = exp ’ (1) 

where A and r are free parameters in the MCMC describing 
the amplitude and timescale of the covariance function, and 
<Ji is the uncorrelated uncertainty for measurement i. The 
log likelihood C for the residuals r is calculated as: 

log£= -ir^E’^r- i log | E | log(27r). (2) 

The Gaussian process modelling provides more realistic un¬ 
certainty estimates. The A and r factors are presented in 
Table |6] 

3.1.2 Results 

The final masses and radii are Mi = 0.244^q qq 3 Mq, Ri = 

0.2611^09 Re, M 2 = 0.179lHo? Re, R 2 = o'.216tH?I Re- 

The full set of system parameters are presented in Table |6l 
The results presented are the mode of the posterior proba¬ 
bility distributions, with the uncertainties representing the 
68 % confidence region. 

The posterior probability distributions for the free pa¬ 
rameters in our global fit are shown in Figure |6l For each 
pair of parameters, we calculate the Pearson R correlation 
coefficient to check for degeneracies. We note a signihcant 
degeneracy between the light ratio L 2 /T 1 and the radius ra¬ 
tio R 2 /R 1 in our results. For systems with partial eclipses. 
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Figure 5. A 35 day sample segment of the HATSouth light curve residuals to the best fit model from Section |3.II The data points are 
plotted in grey, 50 point-width boxcar smoothed function in black, and the cosine filtering model by the red line. The times of transits 
and eclipses, masked out in the data, are marked by vertical red dashed and dotted lines, respectively. The corrected data are plotted 
on the bottom panel. 


it is difficult to constrain both parameters simultaneously 
purely from the transit fight curves without a prior on the 
light ratio. We also note a correlation between (i?i + R 2 )/a 
and inclination i, since both parameters together determine 
the transit shape. The time parameters To and e cos ui are 
correlated, since the phase of the secondary eclipse, deter¬ 
mined by ecosw, is influenced by the ephemeris of the pri¬ 
mary transit Tq. We examined the posterior probability dis¬ 
tributions for the mass and radius of both stars, and found 
no correlation between the derived mass vs radius distribu¬ 
tion. 


3.1.3 Testing the robustness of the results 

To test the dependence of our results against various as¬ 
sumptions and datasets, we re-ran the full MCMC anal¬ 
ysis with different subsets of data and initial parameters. 
For each test, we describe the changes to the free param¬ 
eters, assumptions, and datasets used, and report the ma¬ 
jor differences of the derived results to those reported in 
Table [6l While these tests provide a useful guide for gaug¬ 
ing the robustness of our results, the uncertainty estimates 
we report are likely still optimistic, as they are largely de¬ 
pendent on the per-point measurement uncertainties. Only 
repeated velocity measurements and photometric follow-up 
observations will provide the most robust stellar properties 
for HATS551-027. 

• We set limb darkening coefficients as free parameters, 
but found that we can only place very weak constraints 
on the coefficients from the fitting. In R band, component 
A has best fit limb darkening coefficients of ui = O.SIq'j, 
U 2 = 0 . 31 q; 4 , component B of ui = O.tIqJ, U 2 = O.Ojto', 4 , 
and in I band, component A of Mi = M 2 = O.ljlQ'g, 

component B of Mi = O.TIq;®, M 2 = O.iIq;^. We found the 
limb darkening parameters were largely orthogonal to the 


other parameters. The resulting best fit parameter values 
did not differ significantly between the limb darkening free 
and fixed analyses, and the radii and mass uncertainties were 
not noticeably inflated by freeing the limb darkening param¬ 
eters. 

• To test the effect of the single epoch of PFS veloc¬ 
ities on the resulting derived masses, we refit the obser¬ 
vations without the PFS velocities. The resulting orbit 
semi-amplitudes were = 42.4jjQ'2 kms“^ and K 2 = 
58.01q;| kms“^, yielding masses of Mi = 0.249jlQ'QQ5 M© 
and M 2 — 0.18llg QQ 2 Mq. The derived properties are con¬ 
sistent with the global fit results to within errors. 

• To test if the Ha derived velocities, which reflect veloc¬ 
ities of the chromosphere, are representative of the photo- 
spheric velocities, we fit the system using only WiFeS CCF 
derived velocities. We derived velocity semi-amplitudes of 
Ki = 44.lj;5 Q kms“^ and K 2 = 53.212 2 kms“^, yielding 
masses of Mi = 0.2llQ g 2 ^0 ^md M 2 = 0 . 181 q'q 4 Mq. In 
comparison, using only the WiFeS Ha derived velocities, we 
derive semi-amplitudes of Ki = 43 . 3 ll 4 kms“^ and K 2 = 
59.7llg kms“^, yielding masses of Mi = 0 . 261 q'q 2 

M 2 = 0.1910 04 Mq. The discrepancy between the K 2 val¬ 
ues is likely due to the difficulties in resolving and fitting the 
secondary CCF and Ha component using the lower resolu¬ 
tion spectra. The derived semi-amplitudes and masses are 
consistent to the global fit results to within 2a. 

• To test the dependence of the results on the four 
FEROS measurements and associated uncertainties, we fit 
the dataset with only the FEROS velocities and derive ve¬ 
locity semi-amplitudes of Ki — 42 . 41 q '3 kms“^ and K 2 = 
57.8113 kms“^, yielding masses of Mi = 0.247llgj!,g Mq 
and M 2 = 0.182llgg5 Mq. The masses determined are 
consistent, but have three times larger uncertainties from 
the global fit. To check the dependence of our derived 
masses against the measurement precision, we inflate the 
FEROS uncertainties by 2x, and re-fit using only the in- 
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Figure 6. Posterior correlation diagrams for the free parameters in 
of 0.5 < R < 0.8 are coloured in blue, R > 0.8 coloured in red. 

flated FEROS velocities, deriving semi-amplitudes of Ki = 
42 . 4^0 5 and K2 = kms“’^, yielding masses 

of Ml = 0 . 248 lg;Qig Mq and M2 = O.ISIIq ggg Mq. Inflating 
the velocity measurement uncertainty by two times leads to 
a ~ 2 times inflation of the derived mass uncertainties. 

• In Section 13.1.11 we tried to account for spot-induced 
uncertainties by 1) Monte Carlo of out-of-transit baseline 
flux and 2) addition of Gaussian processes modelling. To 
test their effects, we removed the two treatments and re-fit 
the data. The derived radius ratio is R 2 /R 1 ~ 0.82^004, 
light ratio L2/L1 = O.SSlo'gg, resulting stellar radii of 
Ri = 0.260tg;^^^ Rq and R2 = 0.214t°;°°® Rq. Note the 
significant reduction in the estimated uncertainties (~ 60%) 
without accounting for spot activity. 


the global fit. Pairs of parameters with Pearson correlation coefficient 


• To test the dependence of the results on the system¬ 
atic errors associated with each FTS secondary eclipse ob¬ 
servation, and the polynomial normalisation applied to the 
light curves, we fit the observations using each FTS sec¬ 
ondary eclipse observation independently. Using the FTS 
2012-12-12 light curves for the secondary eclipse, we derive 
a radius ratio of R 2 /R 1 = 0 . 84 lg;g®, fight ratio I/2/T1 = 
0 . 55 tQ gg, resulting stellar radii of Ri = 0 . 259 tQ gg® R© 
and R2 = O.22t^;g4R0. Using the FTS 2013 - 03-20 light 
curves for the secondary eclipse, we derive a radius ratio of 
R 2 /R 1 = 0 . 8 llg;Q 7 , light ratio L2/L1 = 0.521°;)’®, resulting 
stellar radii of Ri = 0 . 262 lg;^J^ Rq and R2 = 0.22lS;g2 R@- 
In each case, the derived values are consistent, with the 
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Table 6. HATS551-027 System Modelling Parameters 


Parameter 

Value Source 

Modelled system parameters 
Period (Days) 

4.077017 (ti)' 

GM*’ 

To (BJD) 

2456374.0171 (tl 

) GM 

(i?l + R 2 )/ O' 

0.0590 hi 

) GM 

R 2 IR 1 

0.82 hi 

) GM 

in 

87.99 (^3 

) GM 

L 2 /L 1 (6114-6574 A) 

0.56 (t®^ 

) GM-I-PFS 

L 2 IL 1 (i?-band) 

0.53 (t^^ 

) WiFeS 

L 2 IL /1 (7-band) 

0.50 (^3 

) WiFeS 

e cos a; 

0.00013 (tg 

) GM 

e sin a; 

-0.002 (tj 

) GM 

Ki (kms“^) 

42.4 (t^^ 

) GM 

K 2 (kms“^) 

57.6 hi 

) GM 

7?-band star A ui 

0.3762 Fixed'^ 

72-band star A U 2 

0.4574 Fixed‘S 

72-band star B ui 

0.5392 Fixed'^ 

72-band star B U 2 

0.3485 Fixed‘S 

2 -band star A ui 

0.1539 Fixed= 

2 -band star A U 2 

0.6330 Fixed'^ 

2 -band star B ui 

0.3124 Fixed'^ 

2 -band star B U 2 

0.5509 Fixed= 

Gaussian process modelling parameters 


A (FTS 20121212) 

0.0005 (tl 

) GM 

T (FTS 20121212) (days) 

0.0013 (^3 

) GM 

A (Swope 20130226) 

0.0022 (t® 

) GM 

T (Swope 20130226) (days) 

0.0037 (^2 

) GM 

A (FTS 20130320) 

0.0004 (t^ 

) GM 

T (FTS 20130320) (days) 

0.0019 (tj 

) GM 

Derived parameters 

Ml (Mq) 

0.244 hf 

) GM 

Ri (Re) 

0.261 (tg 

) GM 

M 2 (Mq) 

0.179 (ti 

) GM 

R 2 (Rq) 

0.218 (15^3^ 

) GM 


“ Uncertainties in parenthesis for the last significant figure 

^ GM - Global modelling of light c urves and radi al velocities 

^ Fixed to values interpolated from lcia^ 3 ll2000l^ using the PHOENIX models. 


uncertainties for the radius and light ratio parameters in¬ 
creased by ~ 30 %. 

• In the Swope light curve, four measurements at phase 
~ 0.15 were consistently below the out-of-transit median. 
We masked out the outlying measurements and performed 
the refit, deriving a radius ratio of R 2 /Ri ~ 0.841q q|, 
light ratio L 2 IL 1 = O.SSIqqy, resulting stellar radii of 
Ri = 0.2611;^;22? i ?0 and R 2 = 0.22l“;°J i?®. The differ¬ 
ences are negligible compared to the global fit. 

• We note that the follow-up photometric observations 
were all obtained in the i band. Single-band follow-up ob¬ 


servations can bias the derived radius ratio parameters if 
the light ratio is band-dependent. For comparison, we fit 
the observations using the r-band HATSouth observations 
only, without the photometric follow-up datasets, we derive 

To = 2456374.017jlo oo 2 (consistent, but 10 times worse in 

precision), {R\+R 2 )/a = 0.0601q qq 2 (consistent but 4 times 
worse in precision), i? 2 / 7 ?r = 0.83jlQ Q 5 (consistent and equal 
precision), light ratio of L 2 /L 1 = 0.521qq 7 (consistent and 
equal precision), and component radii of Ri — 0.267lg Q 52 
and R 2 = 0.22lg qJ (a larger R\, but still consistent within 
Icr). 
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• The light ratio L 2 /L 1 in our global fit is tightly con¬ 
strained by a Gaussian prior of mean 0.51 and standard 
deviation 0.04, determined by the analyses of the PFS spec¬ 
tra. To understand the effect of this prior on our result¬ 
ing parameters, we modify to the Gaussian prior to have a 
mean 0.51 and standard deviation 0.08 (2x broader). The 
resulting fit has a R 2 /R 1 = O.SSjlQ Qg, light ratio I/ 2 /T 1 = 
0.641q J 4 , resulting stellar radii of Ri = 0.25ljlQ Qjg Rq and 
R 2 = 0.231^2 Rq. The results remain consistent with Icr, 
with the radius uncertainties increased by a factor of 1.5. 


3.2 Spectral Classification 

3.2.1 Effective Temperature 


We use the WiFeS low resolution spectrum, which has the 
largest spectral coverage and highest S/N, to constrain the 
effective temperature of HATS551-027. Whilst the spectrum 
of HATS551-027 is a composite, at low resolution we can¬ 
not resolve the two binary components by their velocity 
differences. We first treat the spectral classification as if 
HATS551-027 is a single star. 

The flux calibrated WiFeS R — 3000 observation 
is matched to sy ntheti c spectra from BT- Settl models 
llAllard et al.|[201 with lAsplund et al.l (l2009l l abundances. 
We fix log g to 5.0 (derived from Section[3TT2]to be 4.99 and 
5.01 for star s A and B, and con sistent with stellar isochrone 
models, e.g. iDotter et al.ll200^ and [M/H] to 0.0, and min¬ 
imise the fit by varying Teff at steps of 100 K, interpolat¬ 
ing the Teff-X^ relationship with a B-spline. 

First, to understand the uncertainties in our temper¬ 
ature classification, we obtained WiFeS spectra of five Id- 
dwarfs with temperatures measured from interferometry 
observations. The WiFeS BT-Settl model fit temperatures 
and the interferometry temperatures are shown on Ta¬ 
ble [T] We find an error of 85 K encompasses 68% of the 
deviations between the WiFeS measurements and litera¬ 
ture values, weighted by the literature uncertainties. We 
find an error of 15 0 K when we use models employing the 
I Caff an et al.l (l201lf ) abundances. We also tested the stan¬ 
dard stars against mod els with metallic i ties o f [M/H] = -0.5, 
+0.3, and +0.5 fusing lAsplund et al.l ll2009l i abundances). 
We found a mean error of lOOK, 85K, and 95K for each 
of the metallicities respectively. The temperatures derived 
using the different metallicities differed systematically from 
that of the temperatures derived from the [M/H]=0.0 mod¬ 
els by -40K, +20K, and +60K respectively. 

For HATS551-027, we find a best fit binary effective 
temperature of Tefj = 3113 K. We discuss metallicity later 
in the section. The WiFeS spectrum and the best fit model 
are plotted in Figure [T] 

For comparison, we also plot in Figure [3 the WiFeS 
spectra of GJ 551 (M5.5V 309 8 +56 K, Demory et al.ll2009l l 
and GJ 699 (M4V 3224+lOK. lBovaiian et al.ll2012l 'l'~two M- 


dwarfs with similar temperatures and spectra as HATS551- 
027. A visual inspection of some temperature sensitive fea¬ 
tures, such as the depth of the TiO absorption bands, the Na 
and CaOH lines, confirms that HATS551-027 has a spectral 
type between the GJ 551 and GJ 699. 

To estimate the temperatures of the two stellar compo¬ 
nents, we simultaneously fitted for the temperatures of both 
stellar components using a grid of simulated composite spec¬ 


Table 7. WiFeS M-dwarf Spectral Classifications 


Object 

Literature T^ff (K) 

WiFeS Tbs (K)“ 

Ref 

GJ 191 

3570 ± 156 

3787 

1 

GJ 551 

3098 ± 56 

3000 

2 

GJ 699 

3224 ± 10 

3305 

3 

GJ 880 

3713 + 11 

3671 

3 

GJ 887 

3797 ± 45 

3724 

2 


References l.lSegransan et al.l ll2n0.‘4l 2.lDemorv et al.l (120091') 3. 

iBovaiian et alj 1201^ -^ ^ 

° Matched aga i nst B T-Settl spectral models using 
lAsplund et al.l ||2009^ abundances, with [M/H]=0.0 


tra. We generated a 2D grid of composite spectra from the 
BT-Settl models at steps of 100 K for both simulated stel¬ 
lar components, with fluxes scaled using the best-fit radius 
ratio, and wavelength shifted in velocity according to the Ke- 
plerian orbit solution (both derived in Section [3.HI . We then 
calculate the x^ of fit to the WiFeS spectrum and WiFeS 
CCF-derived light ratios in the R7000 and 17000 spectral 
bands fSection l2.3.1ll . The best fit component temperatures 
were Teff,i = Teg ,2 = 3100, consistent with the binary effec¬ 
tive temperature. 

We can also approximate the temperature of the two 
stellar components by assuming they follow the Stefan- 
Boltzmann law. We adopt a binary effective temperature 
of 3114 ± 85 K, light ratio of 0.50 /)q q 4 , and the stellar radii 
from Table |6l The derived temperature ratio is 0.94 ± 0.04, 
with the temperature of the primary stellar component as 
Teff.i = 3190 ± 100 K and the secondary component as 
Tbh ,2 = 2990 ± 110 K. We adopt these temperature ratio 
derived temperatures for the further analysis. 


3.2.2 Metallicity 


Accurate metallicities for M-dwarfs are difficult to mea¬ 
sure from optical spectra, since the molecular absorption 
features distort the continuum significantly, and mask out 
the common metal lines. We first use the ((TiO/CaH index 
to estimate the metallicity of HATS551-027. We measure 
the band indices for T i05, GaH2, and C aH3, according to 
the regions defined in iReid et al.l lll995l l. to be 0.36, 0.39, 
and 0.68 respectively. We calculate a series of Q values us¬ 
ing relations between C aH2+CaH3 a n d the [Ti05lzFi in¬ 
dex fr o m calibrations in Leoine et al.l (l2007l ): TDhital et al.l 
ll2012ll : lLepine et al.l (l2013li . all yielding ^ = 1.0, the average 
abundance for disk stars. The derived metallicity estimate 
for HATS551-027 is [M/H] = 0.0 ±0.2 using the calibrations 
and u ncertainties between and [M/H] from iMann et al.l 
ll2013ll . For a consistency check, we also use t he metallic¬ 
ity, J — K and V — K colour relationships from iMann et all 
ll2013ll to calculate a photometric colour-derived metallicity 
of [M/H] = +0.1 ± 0.2, consistent with the band indices 
result. However, we note that infrared spectroscopy will be 
necessary to provide a definitive measurement for the metal¬ 
licity of HATS551-027. 
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Figure 7. WiFeS R = 3000 spectra of HATS551-027 and two standard M-dwarfs of similar a spectral type, GJ 551 and GJ 699. The 
observations are plotted in black, best fit BT-Settl spectral model to HATS551-027 in red. 


3.2.3 Activity and Ha Emission 


the two stellar components. 


The low resolution WiFeS spectrum of HATS551-027 shows 
all the Balmer lines, and the Calcium H & K lines, in 
emission. The strength of Ha emission is commonly inter¬ 
preted as an indicator for stellar activity, and potentially 
corre l ated to the inflated r adii of M-dwarfs llLopez-Morales! 
I 2 OO 7 I; IStassun et al.| I 2 OI 2 I) . We use the PFS spectra to 
measure the equivalent width of the Ha emission for both 
stellar components of the binary, fitting Gaussian profiles 
to each feature. We calculate the equivalent width by 
assuming a polynomial fit to the continuum flux, scaled to 
the relative light contribution from each stellar component 
(L 2 /I /1 = O.Sbj^Q Qg Section [3.111 . The measured equivalent 
widths are 2.8 ± 0.2l and 3.6 ± 0.41 for HATS551-027A 
and B, respectively. To calculate the Lhc/LboI ratio for 
each stellar component, we measure a LboI for each stellar 
component by integrating under the respective best fit 
BT-Settl model spectra. The Lhoc/LboI ratio is calculated 
by comparing the flux of the Gaussian fit to the Ha lines to 
the integrated fluxes of the BT-Settl models. The largest 
error contribution comes from the TboI estimate, stemming 
from the uncertainty in the atmospheric parameters. To 
estimate the uncertainty, we calculate a series of Lhoc/LboI 
values for each component whilst adopting a grid of spectral 
models, with Tefr.i from 3000 to 3200 K, and TeB ,2 from 
2900 to 3100 K. The standard deviation scatter in the 
resulting Lhci/LboI is adopted as the uncertainty. The 
derived Lhc/LboI are —4.15 ± 0.03 and —4.08 ± 0.05 for 


The set of stellar photometric and spectroscopic param¬ 
eters for HATS551-027 are presented in Table [51 


4 DISCUSSION 


Following CM Draconis and KOI-126, HATS551-027 is the 
third well characterised double lined eclipsing binary with 
both components in the fully convective regime. Figure [8] 
shows the masses, radii, and temperatures of the HATS551- 
027 system with respect to other well characterised double- 
lined eclipsing binaries in the same mass regime. Figure [8] 
also compares the prope rties of HATS551-0 27 to isochrones 
from the Dartmouth (|Dotte^^£_^ 200^)j_ Bax^e_^£_^ 


PARSEC models 


Chen et al 


,_ _ ■_ ressan et al.l I2012I-'_ 

201 ^, and Yonsei-Yale lYY. IyI et al.l I2001I : l^ada et al 
2 OI 3 I 1 models. 


The radius and estimated temperatures of these bi- 
naries are largely d iscrepant to the Dartmouth, YY, and 
iBaraffe et al.l ( 1998l l models, but generally agree with the 
PARSEC model predictions to ~ la. Compared to the 
Dartmouth 5 Gyr isochrones, the observed radii for the 
HATS551-027 system are larger by 2.6±3.5% and 9.1±4.8% 
for components A and B respectively; the temperatures are 
cooler by 3.8 ± 3.6% and 7.9 ± 3.7% respectively. The mea¬ 
sured system properties are more consistent to the PAS- 
CAR 5 Gyr isochrones, the observed radius of component 
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Table 8. HATS551-027 Stellar Parameters 


Parameter 

Value 

Source 

RA (J2000 HH:MM:SS.SS) 

05:44:57.92 

2MASS 

DEC (J2000 DD:MM:SS.SS) 

-24:56:09.6 

2MASS 

Proper Motion (mas/yr RA DEC) 

229 -25 

Salim & Gould f2003^ 

Proper Motion Error Ellipse 
(mas/yr mas/yr Position Angle) 

5 5 104 

Salim & Gould ('20031 

Photometric parameters 

B (mag) 

17.3 ±0.1 

APASS 

a' (mag) 

16.3 

APASS 

V (mag) 

15.54 ±0.05 

APASS 

r' (mag) 

15.0 

APASS 

i' (mag) 

13.221 ± 0.05 

APASS 

J (mag) 

10.75 ± 0.02 

2MASS 

H (mag) 

10.13 ± 0.02 

2MASS 

K (mag) 

9.85 ± 0.02 

2MASS 

Spectroscopic parameters 

Teff.l (K) 

3190 ± 100 

WiFeS 

Teff,2 (K) 

2990 ± 110 

WiFeS 

[M/H] 

0.0 ±0.2 

WiFeS 

log Tjja.l/Tbol.l 

-4.15 ± 0.03 

PFS 

log Trfa, 2 /Tbol ,2 

-4.08 ± 0.05 

PFS 


A is smaller than the models by 0.6 ± 3.6%, and compo¬ 
nent B larger than the model by 4.2 ±4.9%, both consistent 
within errors. The temperatures for both components are 
also consistent to within errors, with A and B hotter than 
the models by 3.1±5.0% and 2.6±5.6% respectively. The un¬ 
certainties are derived via a Monte Carlo exercise, drawing 
the observed sample from Gaussian distributions based on 
the measurement errors in mass, radius, and metallicity. For 
trials with sub-solar metallicities, we interpolate between 
the [M/H] = —0.5 and 0.0 isochrones to derive the model 
radius, for super-solar metallicites, we interpolate between 
the [M/H] = 0.0 and ±0.5 isochrones. 

For tidally locked binaries, a shorter orbital period, and 
thereby rotation period, may lead to more powerful internal 
dynamos. In turn, the higher stellar activity may result in 
a greater discrepancy between the observe d properties of 
these binaries and their model predictions. [Fleming et al.l 
lll989ll established a rotation - X-ray activity correl ation 
for main sequence single stars. iLopez-Moraled (l2007h and 
IStassun et al.l 1I2OI2I ') demonstrated a correlation between X- 
ray activity and the radius, discrepancy of low mass 
binaries. 


HATS551-027 has a significantly longer period (4.1 
days) than both CM Draconis and KOI-126 BC (< 2 days). 
However, there appears to be no direct correlation between 
the model discrepancy and orbital period. In the HATS551- 
027 system, only HATS551-027B is significantly inflated. In 
the short period (1.3 days) CM Draconis binary, both stellar 
components a re discrepant from th e models in radius and 
temp erature (iTorres fc Riba^ I 2 OO 2 I: iMacDo nald fc Mullani 
\201^ : ISnada et al. 20131 : Feiden &: Chabov^ 2014ll . Com¬ 
pared to the Dartmouth 5 Gyr isochrones, t he radii of the 
CM D raconis stars (adopting parameters from iMorales et ^ 
I 2 OO 9 I I are inflated by 5.5±1.0% and 6.1±1.1%, and temper¬ 
atures cooler by 8.0 ±2.3% and 7.6±2.3% for components A 


and B respectively. However, standard isochrone models can 
reproduce the radi i of the KOI-126 BC system very well (e.g. 
I Feiden et al.l I 2 OIII : ISnada fc Demara~u3 2012l 'l , even though 
it is a similar close-in binary system with period of 1.8 days. 
Compared to the Dartmouth models, KOI-126 B is negli¬ 
gibly smaller by 0.2 ± 1.2%, and C is negligibly larger by 
0.5 ± 1.3%. 


A lack of correlation remains when we also consider 
well characterised double-line eclipsing binaries with one 
low mass component. LSPM J1112±7626 (lirwin et al.ll20lih 
has a long period of ~ 41 days, with both components of 
the binary significantly inflated in radius by ~ 5%, and 
cooler in temperature by ~ 10% compared to standard mod- 
els. Kepler-16, ano ther 41-day period unequal mass binary 
(iDovle et al.ll201ll ') is also inflated in comparison to mod¬ 
els by_~_7%_Twe_^opt the dynamical masses measured 
by Bendei_et_^ J 2 OH ). This is consistent with the findings 
of iFeiden fc Chabovei (I2OI2I '). which showed a lack of pe- 
riod correla t ion fo r detached double lined eclipsing binaries. 

I Zhou et al.l (l2014l ~) also examined a larger sample of double 
lined and single lined binaries in < 0.3 Mq regime, finding 
no correlations with period. 


We also examine Hq, as a direct activity indicator, 
for correlations against the radius and Tefj model discrep¬ 
ancies. The stellar components of HATS551-027 show Ho 
emissions and spot modulation, signatures of significant stel¬ 
lar activity. Both components of HATS551-027 exhibit dis¬ 
crepancies in the Teff with respect to the Dartmouth mod¬ 
els. CM Draconis sy stem is also known to be active from 
it spo t modulations jMor ales et al]l2009ll and X-ray emis¬ 
sion (iLopez-MoralesI l2007l l. Although we lack a directly 
measured Ha emission of CM Draconis, we can estimate 
the Ha flux by converting the X-ray fluxes of CM Draco¬ 
nis reported by iLopez-Morale^ (l2007ll to log Lira/Lboi us¬ 
ing the empirical relationship established for field M-dwarfs 
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Figure 8. The measured masses, radii, and temperatures of the well characterised dou ble-line eclipsing binaries with respe ct to the 
model isochrones. Sola r metallicity, 5 Gyr mo dels from Dartmouth llDotter et a |200St) are plott ed in black, Yonse i-Yale llYi et alj 

in green, and fearaffe et al.lll99^ in blue. Stellar 


I 2 QOII : [s^pada et al.ll2015ll in magenta, PARS EC llBress an et a lj|2012l ; fchen et aPl^^H 


properties for CM Draconis are taken froml^orales et alj l ^Osh. K OI-126 BC from Carter et_ajj _j20n h the ra dii of Kepler-16 B from 
I Doyle etajj jgOjd^ and dynamical masses from lBend^ et alj ||2012^ . LSPM1112+7626 B from llrwin et alj ll201lh . and WTS 19g-4-02069 
B from [NefeeT*^ ll2013h . 


in IStassun et al.l (l2012li . The approximate Ha fluxes are 
log L/fa/Tboi = —S.Sjij g for CM Draconis A, and —3.7jij g 
for B, similar to the Ha activity we mea sure for HATS551- 
027. WTS 19g-4-02069 (iNefs et al.ll2013h was also reported 
to exhibit strong Ha emission and Teg-model discrepancy, 
but we could not estimate log Lira/Tboi value from the re¬ 
ports. One exception to this trend is LSPM J1112-I-7626B, 
which has no reported Ha emission Jlrwin et al ][2oi3) , but 
an estimated temperature discrepant with the Dartmouth 
at 6%. 

We also measure the log I/rrc/I/boi for the single stars 
GJ 191, GJ 551, and GJ 699 from the WiFeS spectra that 
we took for spectral classifications ISection 13. 2. 111 . GJ 551 
exhibits a distinct Ha emission feature, but we can only 
provide Icr upper limits for the Ha strength of GJ 191 and 
GJ 699 (Tabled. 

Figure [9] plots the log L_ffa/Tboi values of HATS551- 
027 and the single stars against the temperature and ra¬ 
dius discrepancy of the respective stars. For reference, we 
also p lot the empirical correlations found by IStassun et al] 
ll2012h . Whilst the small sample size (3) prevents any 
meaningful statistical interpretation, we note that the T^a- 
log L/fa/Tboi relationship appears to be consistent to this 


Table 9. Ho emission strength measured from WiFeS spectra 


Object 

log Lho. /Tbol 

GJ191 

< -4.69 

GJ551 

-4.13 ±0.18 

GJ699 

< -4.63 

HATS551-027A 

-4.15 ±0.02 

HATS551-027B 

-4.03 ±0.05 


analysis for the VLMSs. The two stars without apparent Ha 
emission, GJ 191 and GJ 699, have effective temperatures 
hotter than the three stars with Ha emission features, GJ 
551, HATS551-027 A and B. This is potentially indicative 
that the M-dwarfs with cooler photospheric temperatures 
have hotter chromospheric temperatures. 

Given the smaller model-observation discrepancy in ra¬ 
dius, we cannot confirm a radius-Ha activity relationship 
from this sample. More double-lined binaries with Ha or 
X-ray measurements need to be made to confirm these rela¬ 
tionships. In addition, Ha activity cannot be the sole cause 
of the radius excess, given that HATS551-027A and B ex- 
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Figure 9. The Ha emission strength (log LifaZ-^^bol) for th® five 
VLMSs (TableO are plotted with respect to their radius and tem¬ 
perature deviation from the 5 Gyr solar metallicity Dartmouth 
isochrone. The measurements for HATS551-027 are plotted in red, 
with the two stellar compon ents labelled. We also plot the empiri¬ 
cal relationships derived bv IStassun et al] ll2012fl (S12 in the plot 
labels) using M-dwarfs (< 0.6 Mq) with Ha emission strength 
derived from X-ray fluxes for comparison. 

hibit similar Ha activi ty, but only B ex hibit a radius ex¬ 
cess to the isochrones. lirwin et alJ (l201ll ') also note a lack 
of detectable Ha emission features for LSPM J1112+7626, 
despite the inflated radii for both components of the binary. 

The HATS551-027 system is also unique as the only 
very low mass M-M binary well positioned for southern 
hemisphere facilities. This should enable a host of follow-up 
opportunities to better characterise the properties of this 
benchmark system. 
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